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Spin-currents and non-abelian gauge potentials in electronic systems can be treated by spin- 
current-density functional theory, whose main input is the exchange-correlation (xc) energy ex- 
pressed as a functional of spin-currents. Constructing a functional of spin currents that is invariant 
under U(l)xSU(2) transformations is a long-standing challenge. We solve the problem by express- 
ing the energy as a functional of a new variable we call "invariant vorticity" . As an illustration we 
construct the xc energy functional for a two-dimensional electron gas with linear spin-orbit coupling 
and show that it is proportional to the fourth power of the spin current. 

PACS numbers: 31.15.E-, 31.15.eg, 72.25.-b 



I. INTRODUCTION 

In the last few decades, density functional theory 
(DFT)i has grown to be a widely used method for 
studying the ground-state properties of interacting many- 
electron systems and the range of its applications has 
been expanding. An important generalization of DFT, 
which takes into account the orbital effects of an exter- 
nal magnetic field, is the non-relativistic current density 
functional theory (CDFT), which was formulated by Vi- 
gnale and Rasolt (VR)2~— in the late eighties. The basic 
variable of that theory (in addition to the usual parti- 
cle and spin densities) is the paramagnetic current den- 
sity j p (r) , which has the following advantages upon the 
"physical current density" j(r): (i) it has no explicit de- 
pendence on the external vector potential - a property 
it shares with previous relativistic formulations^, and 
(ii) it does not vanish in the limit of uniform density 
and magnetic field, which is vital to the construction of 
a local density approximation (LDA). Unfortunately j p 
is not a gauge-invariant variable. VR ensured the gauge 
invariance of the exchange-correlation (xc) energy by ex- 
pressing the latter in terms of the vorticity 

v:=Vxiz, (1) 

n 

which is gauge invariant. While this choice is not unique 
(any gauge-invariant field which is in one-to-one corre- 
spondence to the vorticity would be in principle admis- 
sible), it was shown to be the most natural one, in the 
sense of leading to a local description of the effect of the 
magnetic field in the quasi-homogeneous limit. 

In extending their theory to spin polarized systems^ 
VR proved the analogue of the Hohcnbcrg-Kohn theo- 
rem for paramagnetic spin-currents, and observed that 
the xc energy E xc - now a functional of paramagnetic 



spin currents - should be invariant under local U(l) x 
SU(2) gauge transformations. These are transformations 
in which the wave function is acted upon by the position- 
dependent operator exp[zA(r)], where A(r) is a matrix 
in two-dimensional spin space. However, the behavior of 
the paramagnetic spin currents under these transforma- 
tions is quite complicated and VR were unable to identify 
a variable analogous to the vorticity, such that an E xc 
functional of this variable would be manifestly gauge in- 
variant. 

Recently, Bencheikhi and Gorling^ have further ex- 
tended spin-CDFT to incorporate spin-orbit (SO) inter- 
actions. These interactions are naturally described in 
terms of SU(2) gauge potentials, i.e., vector potentials 
that arc matrices in spin space. 

In Ref. [IH it was noticed that in many realistic situa- 
tions the SU(2) vector potential produced by the SO cou- 
pling has a nonvanishing covariant curl, which is equiva- 
lent to the presence of a non-Abelian magnetic field. This 
effective magnetic field generates "diamagnetic" spin cur- 
rents in the ground state of molecules and solids with 
non-negligible SO interaction. In spite of the fact that 
the ground state spin current is hard to detect experimen- 
tally, its inclusion into the set of basic variables of DFT is 
important to correctly account for the corresponding con- 
tribution to the energetics. SU(2) gauge potentials also 
appear in effective hamiltonians for spin-transfer torque 
systems^, and in pscudospin-orbit coupled systems like 
graphenei^. 

Bcnchcikh derived the form of the Kohn-Sham equa- 
tion in terms of effective U(l) x SU(2) gauge poten- 
tials, both expressed as functional derivatives of a gauge- 
invariant E xc . But again, no explicitly SU(2)-invariant 
expression for E xc was provided. While the main chal- 
lenge of CDFT is to construct practical approximations 
for the xc functional, such as the exact-exchange current- 
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and spin-current density functionals recently developed 
in Refs. [8lfTll and [TH, we believe it is of fundamental 
importance to know how to construct manifestly gauge- 
invariant xc functionals of the spin-current densities. 

In this paper we present a solution to this problem. 
First we note that both the external scalar U and vec- 
tor A potentials and the SO terms can be combined into 
a U(l)xSU(2) four- vector gauge potential, collcctivclly 
denoted by A, in such a way that the exact many-body 
Hamiltonian becomes invariant under local U(l)xSU(2) 
gauge transformations. This leads to the next key ob- 
servation that the ground state energy E depends on A 
only through the invariant field strength B[.A], defined as 
the "SU(2)-invariant curl" of A: all these quantities will 
be defined precisely below. B is invariant under SU(2) 
rotations in spin space and transforms as a pseudo- vector 
in the ordinary space, so we can write 

E[A] = E[B[A}} . (2) 

where E is a scalar functional of B. Notice that in the 
familiar spinless case A reduces to the usual vector po- 
tential A and B reduces to the magnetic field B = V x A, 
so the above equation simply says that the energy is a 
functional of the magnetic field. We will show below that 
the quantity M~ 1 Sp, where 3 v denotes the paramag- 
netic particle/spin currents and N is the particle/spin 
density (again, precise definitions will be given below), 
behaves under U(l) x SU(2) gauge transformations pre- 
cisely like the gauge potential A. This suggests a path to 
constructing a gauge- invariant xc energy functional E xc . 
Let us introduce the "SU(2)-invariant vorticity" 

V := B[N~ 1 J p ] , (3) 

which is the "SU(2)-invariant curl" of N~ 1 J' P and will 
be shown to reduce to the VR vorticity dJ) in the spinless 
case. Since N~ 1 3~ p transforms as A it is evident that 
V is invariant under SU(2) rotations in spin space, just 
as B[A] is. Therefore, the gauge invariance of E xc will 
be guaranteed if we write it as a scalar functional of the 
SU(2)-invariant vorticity (hereafter abbreviated as "in- 
variant vorticity"): 

E XC [J P ] = E XC [V] (4) 

(the dependence on particle and spin density is left im- 
plicit - see concluding section). 

The functional dependence of E xc on V is of course 
very complex and nonlocal. The test of usefulness of 
Eq. is whether there are limiting cases in which the 
E xc functional can naturally be cast in this form. We will 
show that this is the case. Namely, in the important limit 
of electron gas subjected to uniform and weak SU(2) po- 
tentials the xc energy can be written as a quadratic local 
functional of the invariant vorticity. Remarkably, in this 
regime the invariant vorticity is a quadratic functional 
of the spin currents themselves, not of their derivatives. 
This opens the way to a genuine LDA, completely differ- 
ent in character from the local approximation of Refs. 



and [H in which the xc energy depended on derivatives of 
the paramagnetic current and/or density. The leading 
contribution to the present LDA functional goes as the 
fourth power of the spin current, and the xc potentials 
are proportional to its third power. 

This paper is organized as follows. In Section [TT] we 
introduce our model and basic variables. In Section IIIII 
we present the gauge invariant spin-current-density func- 
tional theory and the self-consistent Kohn-Sham equa- 
tions. Our local density approximation for E xc is pre- 
sented in Section IIVI In Section [V] and Appendix [A] we 
calculate the exchange-correlation energy of Rashba sys- 
tem, which would provide the necessary input for our 
LDA scheme in the case of isotropic SU(2) vector poten- 
tials. Finally, Section I VII summarizes our main conclu- 
sions. 



II. BASIC DEFINITIONS AND GAUGE 
TRANSFORMATIONS 

We start with the many-body Pauli Hamiltonian^, 
which describes non-relativistic electrons in a static elec- 
tromagnetic field (h=e = m = c=l) 

H = J dr ¥ | in 2 + i [II • (cr x E) + (er x E) • LI] 

-U • erj* + W , (5) 

where ^ = (^1, ipV) is the two-component field operator, 
n = — iV + A is the kinetic momentum, <r is the vector 
of the Pauli matrices, U is the electrostatic potential, A 
is the vector potential, E = — VJ7 and B — V x A are 
the electric and the magnetic fields, respectively. The 
first term on the right hand side of Eq. ([5]) is the kinetic 
energy, the second is the spin-orbit coupling energy, the 
third and the fourth are the potential and the Zeeman 
energies, respectively. W denotes the electron-electron 
interaction, which we assume to be spin- independent and 
gauge- invariant . 

We can write the hamiltonian more compactly by in- 
troducing the U(l)xSU(2) gauge potential = A^r a , 
where f is a four- vector constructed by a 2 x 2 unit ma- 
trix I and Pauli matrices: (I, a x , <7 v ,a z ) Here and in the 
following upper indices refer to charge [a = 0) and spin 
(a = x 1 y, z), and lower indices to time (fi = 0) and space 
(fi = x,y,z) components. A sum over repeated indices 
is implied^. The components A" of the potential are 
defined as follows 

We also introduce the covariant space derivatives T>j = 
dj + iA'jT , where j = x,y,z, dj is the derivative with 
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respect to the corresponding variable. With this notation 
the Hamiltonian ([5]) takes the form 



H = / dr 



W. 



(7) 



One can easily show that this Hamiltonian is invari- 
ant under a local time-independent U(l)x SU(2) gauge 
transformation 14 ' 16 : 



# — > U^>, 
An — > UAM 1 + HfidJA , 



(8) 



with U = exp [iA a (r)r Q ], where A Q (r) is an arbitrary 
time independent function. This transformation corre- 
sponds to the multiplication by a local phase factor sup- 
plemented with a local rotation in the spin space. Ex- 
plicitly the matrix trasformation law (J5J translates to the 
following transformation for the components A^ of the 
gauge potential 



A a 



K af) {A) 



(9) 



where the trace Tr is taken over r, and TZ(A) is a 4 x 4 
matrix, with TZ 00 = 1, TZ 0a = K a ° = 0, and TZ ab = 
R ab (A), R ab (A) being the 3x3 rotation matrix in spin 
space by angle 9 = 2|A| along the direction A. Note that 
U^doU = for static transformations. 

A crucial quantity for our purposes is the "field 
strength" defined as the covariant curl of the gauge po- 
tential^ 



Bl = B?r a = e m V 5 A%r a 



or, more explicitly 

B?[A] = e ijk [djAt - (1 - 5 a , )e abc A)A c k } 



(10) 



(11) 



The "charge" component (a = 0) of the field strength 
coincides with the usual magnetic field, while its "spin" 
part (a = x, y, z) corresponds to a non-Abelian SU(2) 
magnetic field. Here we must emphasize an important 
difference between the familiar Abelian U(l) magnetic 
field B of electrodynamics and its non-Abelian counter- 
part B. While B is gauge-invariant tout-court, our non- 
Abelian field strength B is gauge- covariant, i. e., under 
the gauge transformation it transforms asKi-> UBlft . In 
other words the behavior of each spatial component Bi 
under SU(2) rotations is that of a vector in spin space: 
Bf -> K al3 (A)B«. The lack of full gauge invariance be- 
comes a problem when we attempt to construct a non- 
local, yet gauge-invariant energy functional, for fields at 
different points in space transform differently. For exam- 
ple, the plausible functional 



dr / dr'K(r,r')B?(r)B?(r') 



(12) 



is not SU(2)-invariant, in spite of the fact that all indices 
are duly contracted. 




FIG. 1: Two admissible choices of paths connecting the origin 
to two points r = ri and T2 for the evaluation of the path- 
ordered integral in Eq. (|13[) . The solid path is a straight radial 
line, the dashed path proceeds along the x-axis first, and then 
along the y-axis. 



To avoid this difficulty we must go a step further and 
introduce an SU(2)-invariant field, which behaves like a 
scalar under SU(2) rotations. To this end, we define a link 
operator L(r) ( a 2 x 2 matrix) in the following manner: 



L(r) = Pexp f —i / Ai(x)dx 



(13) 



where the integral is taken along some conventionally 
specified path which connects the origin (0) to the point 
r. These paths should be smooth functions of r: for ex- 
ample a set of lines radiating from the origin, or a set of 
"L-shaped" lines (see figure [T|) would be acceptable. 

The symbol P denotes the "path ordering operator" 
along the chosen path that connects to r. In analogy 
with the time-ordering operator, P orders r-dependent 
matrices in such a way that the ones with the "earlier" 
r (i.e., closer to the origin) act before the ones with the 
"later" r. Geometrically, the operator L(r) of Eq. fTS"]) 
can be interpreted as a parallel transport operator - it de- 
scribes rotation of a spinor under a "parallel transport" 
along the specified path in the space with nonzero SU(2) 
connection A\ (see, for example, Ref . [l8T) . It is not dif- 
ficult to show that under SU(2) transformation the link 
operator transforms as follows: 



L(r) 
£t( r ) 



U{v)L{v) 



Then we define the invariant field 13 as follows: 



Bi(r) = L\r)Bi(v)L(r) 



(14) 



(15) 



Making use of the previously derived transformations, it 
is easy to check that B is indeed invariant under local 
SU(2) transformations. The reason for this property is 
that Eq. (fT5")) is a "closed loop" form. In fact, the link 
operators L(r) and L'(r) in Eq. ([13)) describe a parallel 
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transport from the point to r and back, respectively. It 
is also worth noting that our invariant ficd 3 is closely re- 
lated to the "twisted curvature" entering the non-Abclian 
Stokes theorem^. The ground-state energy functional is 
therefore naturally expressed in terms of scalar combina- 
tions such as BfBf, which are also invariant under global 
spatial rotations^. In particular, the expression (|12[) be- 
comes SU(2)-invariant if Bf(r) is replaced by Bf(r). 
The physical spin-current densities 



SH 



SE[A] 



SAZ(r) «Ug(r)' 



(16) 



are SU(2)-covariant in the sense that jf(r) — > 

lZ a " (A)j^(r). Their /i = components form the par- 
ticle/spin density 



Jo 



(r) = s a (r) = (*tr Q *) , 



(17) 



where s (r) = n(r) is the particle density. The /i ^ 
components are particle (a — 0) and spin (a =/= 0) 
currents 

ifW = jd*^*- (2>«*)W]) 



^(r)+A^(r)„4f(r), 



where A^ aa 



A/" 0a = A/" 1 



and all the other 



components are zero. The paramagnetic spin-currents 
are defined as-^ 



(19) 



and we also define jp ( r ) = s "( r )- These paramagnetic 
currents are not gauge-covariant, but transform under 
gauge transformation in the following manner: 



7". 



-M^ Tr [r^BiU] 



(20) 



More importantly for the present purpose, the combina- 
tion vf = {J\f~ 1 ) a ^ 'jpi, which is the analogue of a veloc- 
ity field, transforms exactly like the gauge potentials Af , 
namely 



K«f>}vf - ^-Tr [t^ diU] \ . (21 



It then follows that the covariant curl of Vi 
Vi = Vfr a = e ijk V vd v%r a , 
or, more explicitly, 

Vf = e ijk djV% - (1 - S at o)eij k e abc VjVk 



(22) 



(23) 



transforms as Vf — > H a °(A)V", i.e. it is a vector in 
spin space, just as the spin density. In Eq. (|22[) we have 
introduced T>„ 



dj+ivfT a . 



The gauge co-variance of Vi (or the spin density for 
that matter) poses, as noted above, a problem in con- 
structing nonlocal functionals that are invariant under 
SU(2) transformation. One way to solve the problem is 
to introduce, in parallel with the introduction of B, the 
SU(2)-invariant vorticity: 



Vi(r) = Ll(r)Vi(r)L v (r) , 
where L v is defined, in analogy to Eq. (|13[) . as 



L v (r) = Pexp —i I Vi(x)dx 



(24) 



(25) 



In exactly the same way we can define the SU(2)- 
invariant spin-density 



f(r) = Lt(r)S(r)L„(r). 



(26) 



where s(r) = s a (r)r a is the ordinary (covariant) spin 
density, defined in Eq. (fl~7|) . Curiously, the expressions 
for the invariant vorticity and spin density, although in- 
variant, contain the non- invariant paramagnetic currents 
through ii(x). We will return to this point in the conclu- 
sions. We now have all the ingredients for constructing 
a SU(2)-invariant xc energy functional. 



III. SPIN-CURRENT-DENSITY FUNCTIONAL 
THEORY 

The ground state energy and spin current distribution 
of an interacting many-body system, subjected to exter- 
nal gauge potentials A'Z, is obtained from the minimiza- 
tion with respect to of the functional 

E GS [j^} = min MH\t/>) = F\j^] + E ext [j^} . 

(27) 

Here F\j"] is the part of the ground-state energy that 
does not depend on the external fields. This func- 
tional is traditionally expanded as F = T s + En + E xc , 
where T s , En and E xc are the non-interacting kinetic, 
Hartree and exchange-correlation contributions, respec- 
tively. E ext \jp ] , on the other hand, is the part of the 
ground-state energy that pertains to the coupling of the 
currents with the external potentials. This functional has 
the same form in both interacting and non-interacting 
systems. Since F + E ext , T s + E ext and En are sepa- 
rately gauge invariant, the exchange-correlation energy 
alone E xc = F — T s — Eh has to be gauge invariant as 
well. Unfortunately, the paramagnetic current variable 
jp^ docs not possess this invariance. This becomes a 
problem when one attempts to construct approximations 
in terms of j^. Clearly it would be desirable to devise 
a new variable that automatically guarantees the gauge- 
invariance of an approximate E xc . From the work of the 
previous section wc clearly see that the covariant vorticity 
fulfills our needs. We thus write E xc [jp^\ = E xc [n,s,V] 
as anticipated in Eq. @. Additional requirements of 
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rotational invariance forces E xc to actually depend on 
scalar combinations of V™ as we will see in an explicit 
example below. 

Given an approximate functional of this form the single 
particle Kohn-Sham (KS) equation can be written as£ 



1 

2 ^ 



$A = £A$A, 



(28) 

Here &\{r) = (0A,t(r), 4>\^(r)) T is a two-component 
spinor with (f>\, a (r) being a single-particle wave function, 
and the KS potentials are defined as 

-4°ks,o = < + «h + A Q xcfi - \a^ (24? + A a x J , 



A a — A a A- A a 

Aa _ Aa , Aa 



A a A 



(29) 



in order to reproduce same paramagnetic spin-currcnt- 
densities as the interacting system. Here the exchange- 
crrclation potentials are 



5E xc [n,s, V] 
Sj a 

■>r> u 



and the Hartree potential is 

SEa[n] 



5n 



(30) 



(31) 



The paramagnetic spin-current densities can be readily 
obtained from KS orbitals: 



(32) 



where sums run over occupied orbitals. 



IV. LOCAL DENSITY APPROXIMATION 

While the theory of the previous section is formally 
exact, it is not yet obvious that a local density approx- 
imation can be naturally formulated in terms of V. To 
show that this is indeed the case we follow the approach 
of Ref. 1, i.e. we start from the form of exchange- 
correlation energy and potentials for an electron gas of 
uniform density subjected to small and uniform SU(2) 
potentials 4f . This system can be considered as a refer- 
ence for LDA calculations. The leading order correction 
to the ground-state energy for weak Af i&M 



E, 



GS 



E GS + £ / 

Egs + £ / 



(33) 



V20L Ah AC 



where A is an interaction-dependent constant and Eq S 
is the ^.-independent part of the ground state energy. 
From Eq. ([IT]) we see that the uniform field is related to 
the uniform SU(2) vector potentials as follows: 

(34) 

Notice that at this level of approximation, there is no 
difference between using the covariant field B or the in- 
variant ones B: because all the fields are evaluated at 
the same point in space the link operators cancel in the 
trace. The situation would be different if we had a non- 
local functional, but even in that case the difference be- 
tween B and B can be ignored as long as one works to 
second order in the field strength. Notice, however, that, 
in writing Eq. (|33[) we have assumed that the energy is an 
analytic function of B: the latter appears quadratically 
as required by rotational invariance. 

The form (|33p of the energy functional, combined the 
above relation between B and A, and with Eq. (|16|) , im- 
plies the existence of a uniform equilibrium spin current, 
which cubically depends on the potentials^: 

j? = A L^A 1 ;,^ A*A\A§ . (35) 

Hence the paramagnetic spin current is given by the ex- 
pression 

f pi = A {AlA'A) - 4?4^) - A^4f , (36) 

which, to leading order in the potentials (i.e., to first 
order in SO interaction), reduces to the following simple 
form 

#< « -A^4f . (37) 

It is clear from this expression that the field strength co- 
incides (at this level of approximation) with the covari- 
ant vorticity. Thus, the ground-state energy functional 
for weak fields is obtained simply by replacing Bf —> Vf 
in Eq. (|33[) and the corresponding xc energy functional 
is 

4_^E£ / drTr{V,(r)V 4 (r)} 



E r 



E„„ 



A, 



drVf(r)Vf(r), 



(38) 



where E® c is the xc energy of the system with zero vor- 
ticity, and Vf is given, in the limit of slowly varying den- 
sities and currents, by the last term of Eq. (f2"3")l. i.e., 



Vf ~ —eijk e abcVjVk (all indices spatial), and X xc = A — Ao 
is the difference between the fully interacting A and the 
noninteracting one. From Eq. (|38[) the xc potential can 
be easily calculated according to Eq. ([3H)l . For example, 
in the absence of spin polarization we get 

A 

,4 



A a = '— (j a n b i b - j a i b j b ) CW) 
jjA \JpiJp jJp 7 JpjJpiJpjJ ' 



While Eq. {3SJ is of the same form as the first CDFT 
functional proposed for weak fields in Refs. [U and H, we 
point out that the quadratic functional dependence of the 
covariant vorticity on paramagnetic current produces a 
completely different and new dependence of A xc on jjr^. 
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FIG. 2: Excess exchange-correlation energy AE XC (filled cir- 
cles) and excess exchange energy AE X (empty rectangles) of a 
two-dimensional electron gas at Wigner-Seitz radius r a — 1 vs 
Rashba spin-orbit coupling constant a = a/vp (vf = y/2/r B 
is the Fermi velocity) in units of E xc = —0.8 and E x = —0.6 
(in the atomic units) respectively. Inset: same as in the main 
figure but in units of E xc a 4 and E x a 4 , respectively for the 
exchange-correlation and exchange energies. 



FIG. 3: \ X c of a two-dimensional electron gas with Rashba 
spin-orbit coupling vs Wigner-Seitz radius r s , for varying 
Rashba coupling constant a — a/vp (up = y/2/r s ). The 
limiting behavior at a = (solid line) is obtained by extrap- 
olation. 



VI. DISCUSSION 



V. CALCULATION OF X xc IN A 
TWO-DIMENSIONAL ELECTRON GAS 

Let us finally outline the calculation of X xc for a uni- 
form two-dimensional electron gas with linear spin-orbit 
coupling (Rashba modcl22). The model in question^ is 
characterized by a gauge potential with components 

A = -A = a . (40) 

where a is the Rashba coupling constant. This implies 
that 

Bf = -2a 2 , (41) 

is the only non-zero component of B, and from our previ- 
ous discussion we expect AE XC = E xc — E® c oc a 4 . Some 
details of the calculations are provided in Appendix 
An interesting observation is that the exchange energy 
alone, while clearly gauge invariant, is a non-analytic 
function of a, being cx a 4 In a at small q*21. So it is essen- 
tial to include the correlation energy in order to restore 
the expected analyticity of the energy functional. This 
we have done with the help of the random phase approxi- 
mation (RPA)i. In Fig.[2]we plot the ratio AE XC /E° C vs 
a as well as the corresponding results for the exchange 
only. It is immediately evident that AE XC scales as a 4 
but AE X does not. In Fig. [3] we plot the constant X xc as 
a function of Wigner-Seitz radius r s . This quantity will 
be the necessary input for LDA calculations within our 
spin-currcnt-dcnsity functional formalism. 



The results of this paper are primarily formal. The 
introduction of the "invariant vorticity" allowed us to 
provide a formal solution to a long-standing problem in 
spin-CDFT, namely, the problem of constructing man- 
ifestly SU(2)-invariant functionals of the paramagnetic 
current densities. An interesting feature of the present 
construction is that the use of the paramagnetic spin cur- 
rent as basic variables is mandatory: without them we 
could not define the link operators L v (r) (see Eq. ([25)) ). 
which arc vital to the definition of the invariant vortic- 
ity. This is at odds with spinless CDFT, in which the 
use of the paramagnetic current density as a basic vari- 
able is, in a sense, optional. That is, while the use of the 
paramagnetic current is strongly indicated for practical 
reasons (e.g., for constructing a local density approxi- 
mation), one might nevertheless chose with work with 
the physical current, provided one is willing to accept 
a much higher degree of nonlocality in the functionals. 
The situation is quite different here. It is precisely the 
nonlocality of the functionals that forces us to introduce 
the link operator, and hence the paramagnetic current. 
However, this complication does not show up within the 
local density approximation. 

Even more striking is the fact that the paramagnetic 
spin current must necessarily enter the construction of a 
nonlocal SU(2)-invariant functional of the spin density. 
In other words, such a functional must be constructed in 
terms of the invariant spin density s(r) of Eq. (|26p . rather 
than the ordinary spin density. The physical reason for 
this is that SU(2) symmetry cannot be maintained with- 
out the introduction of SU(2) vector potentials, which in 
turn modify the paramagnetic spin current densities. A 
truly SU(2)-invariant nonlocal spin DFT must necessar- 
ily be a spin-current DFT. Even if we stop at the "semi- 
local" level, e.g. by doing a gradient expansion in the spin 
density, SU(2) invariance requires that the gradients be 
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replaced by covariant gradients, which are most naturally 
defined in terms of the connection v^. Another possibil- 
ity is to resort to functionals of spin-orbitals, which can 
be SU(2) invariant without explicitly involving spin cur- 
rents. 

In spite of the formal progress, there is a major weak- 
ness in our proposal. The link operators depend, for their 
definition, on an arbitrary choice of linking paths (see 
Fig. [T]). While it is evident that physical results cannot 
depend on the choice of the paths, we cannot exclude the 
possibility that the form of an approximate functional 
might depend on this choice, even though the functional 
itself is SU(2) invariant for a given choice. This question 
clearly deserves further consideration. 

In conclusion, we hope that our results will stimu- 
late further work aimed at constructing novel exchange- 
correlation functionals for spin-orbit coupled systems. 
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energy per particle, due to spin-orbit interaction reads^ 



AE X = - — -2 / dqqv q da; [U a (q,icj) - IL (q,iuj)] , 
(27r) Jo Jo 

(Al) 

and the excess exchange-correlation energy per particle, 
within the random phase approximation is given by 



AE. 



RPA 



S 

(2^ 



dtu In 



1 - v q Ii a (q,iu)) 



q Jo * [1 - v q Ilo(q,iuj) 

(A2) 

In Eqs. (|A1[) and (|A2[) . S is the area of system, Ho(q, iui) 
is the noninteracting density-density response function 
of a two-dimensional electron gas (2DEG) without spin- 
orbit coupling (i.e. Stern function^), and H a (q,ioj), 
the noninteracting density-density response function of 
a 2DEG with Rashba spin-orbit coupling is 
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n Q (q, iu) 



S ^ ILJ 



^-•^r k+q > (A3) 



-k,^ <--k+q,f 

where is the Fermi distribution function and £k,p 
k/2 — [j,ak. Moreover 



•^k.k+q 



1 + fll/ COS 6>k,k+q 



(A4) 



Appendix A: Details for the calculation of A xc 

Exchange-correlation energy is obtainable by means of 
the fluctuation-dissipation theorem. The excess exchange 



is the form factor, with #k,k+q being the angle between 
k and k + q. After changing sum over k in Eq (|A3jl into 
integral and performing sum over v, we find 



II Q (a;, iu) = — 



dyy d<p- 



dfj, cos < 



a cos z 



6 M cos (j) + c p 



(A5) 



Here y = k/kp, x = q/(2kp), u = u/(k-pq) 1 fi^ = 
^F/i/^F = Vl — a 2 +jjLa and a = a/k-p, with kp = \f2im. 
Here we have also introduced the following notations: 



Integration over (f> can be performed analytically and the 
result reads 



xy 2 , 



b ll = y [2x (x — iu) + a (fj,y — a)] , 



x [x — iu) 



+ iiay (x — iu) , 



dfj. = y + \ia , 

= x — iu + nay/x . 



(A6) 
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U a (x,iu) = ^Y1 

A 



dy -Re 



+ ^A 



a (Vt 



(A7) 



M A 



with u M± = ± - 4 ac ^ /(2a). merically, we find spin-orbit correction to the exchange 

AT . t-, tt-^i • /rrni j /Tttti, j and exchange-correlation energies of a 2DEG. 

Now using Eq. (|A7|) m Eqs. (|A1[) and (|A2|). and per- b b 

forming the three remaining integrals over y, w and o nu- 
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